Interaction with a field: a simple integrable 
model with backreaction 

Ainaury Mouchet 

"qq" ■ Laboratoire de Mathematiques et Physique Theorique, 

O ■ Universite Frangois Rabelais de Tours — CNRS (UMR 6083), 

S:^ ! Federation Denis Poisson, 

Pare de Grandmont 37200 Tours, France. 
niouchet@phys.univ-tours.fr 

^ ■ July 9, 2008 

CU] PACS: 46.40. Cd, 46.40.Ff, 42.50.Ct, 03.65.Yz. 

rH , Abstract 

Cw ' The classical model of an oscillator linearly coupled to a string cap- 






tures, for a low price in technique, many general features of more realistic 
models for describing a particle interacting with a field or an atom in a 
electromagnetic cavity. The scattering matrix and the asymptotic in and 
out waves on the string can be computed exactly and the phenomenon of 
^ ' resonant scattering can be introduced in the simplest way. The dissipa- 

CO , tion induced by the coupling of the oscillator to the string can be studied 

OO ' completely. In the case of a d'Alembert string, the backreaction leads to 

an Abraham-Lorentz-Dirac-like equation. In the case of a Klein-Gordon 
string, one can see explicitely how radiation governs the (meta)stability 
f^^ ' of the (quasi)bounded mode. 

o 

O ■ 1 Introduction 



m 






Getting back to the spirit of the xix'^ century — when purely mechanical mod- 
els were ubiquitous even for understanding systems involving electromagnetic 
j^ ■ fields — this paper discusses a simple model of an oscillator coupled to a string 

that presents a host of interesting features in a very accessible way. It will be 
presented in detail in § [5] (see figure [TJ . Even though the required technical 
background is on an upper undergraduate/early graduate level, if one remains 
within a classical {i.e. non-quantum) context, this model encapsulates many 
relevant features of physical interest which can be found in more realistic mod- 
els, for instance when considering the interaction of an atom (the oscillator) 
with light (the string). 

First (§121), when the string is infinite, it allows one to gain insight into (or 
to discover for the first time in academic studies) the scattering of waves by a 
dynamical system. It provides the opportunity to introduce some of the ingre- 
dients of scattering theory: when dealing with a one-dimensional mechanical 



system, the notion of asymptotic states (the so-caUed "in" and "out" states) 
is made as simple as possible and the S matrix takes a particular simple form. 
Unlike the transmission of light in a Fabry-Perot interferometer, the resonant 
scattering appears here explicitly in connection with the familiar resonance phe- 
nomenon of a forced oscillator and we can also understand how the motion of 
the oscillator affects in return the external excitation. 

Second (§ 2]), from the point of view of the oscillator, this model consti- 
tutes the most elementary example of radiation. It shows concretely how an 
interaction not only induces a shift in the natural frequency of the oscillator 
(this can already be seen when only two degrees of freedom are coupled) but 
also that the coupling to a large number of degrees of freedom (the string being 
seen as a large collection of oscillators) induces a friction term for the oscillator, 
although no dissipation exists in the system as a whole (for an electric analo- 
gous phenomenon see § 22.6 of Feynman et al., 1970; Krivine & Lesne, 2003). 
Indeed, as far as I know, there are very few places where this model is dis- 
cussed and always in the specialised literature with a d'Alembert string (SoUfrey 
& Goertzel, 1951; Dekker, 1985). Some of its variants (Stevens, 1961; Ru- 
bin, 1963; Yurke, 1984; Dekker, 1984; Yurke, 1986) are introduced precisely for 
studying dissipation at a quantum level. As we will see, a Klein-Gordon string 
allows to keep one discrete mode (a bounded state) without dissolving it in the 
continuous spectrum of the string and therefore allows to mimic the interaction 
of a field with a stable particle, not just a metastable one. This model is parti- 
cularly relevant to see how backreaction works: For instance, for a d'Alembert 
string, the oscillator is governed by an Abraham-Lorentz-Dirac-like equation. 
Besides, we will be able to illustrate precisely the deep connection between the 
resonance and the poles of the S matrix (a major feature in high energy particle 
physics and in condensed matter physics). 

Before we give some guidelines for further developments in the conclusion 
(§©, we will complete our classical study in §[n]by the detailed diagonalization 
of the Hamiltonian and the discussion of the completeness of the basis of modes 
that are used to described the dynamics. This will be the occasion to sketch the 
finite size effects if one wants to use this model for describing an atom placed 
in a cavity. This study also prepares the ground for the quantization which will 
be proposed in a future paper. 

2 The spring-string model 

2.1 Description of the model 

A very thin homogeneous string of linear mass /io is considered to have only 
transverse displacements in one direction. At equilibrium it forms a straight line 
along the cc-axis with the uniform tension Tq (we will never take into account 
the effects of gravity). The string is coupled to a system of mass M with 
one degree of freedom connected to a fixed support with a massless spring of 
stiffness K (see figure [T]). The coupling is modelised by a second massless 
spring of stiffness k attached to the string at the massless point A a.t x — 0. All 
the vibrations will be considered within the harmonic approximation of small 
amplitudes. We will denote by ^{x,t) the transverse displacement at time t of 
the string element located at x at equilibrium. The displacement of M with 




Figure 1: Model of a harmonic oscillator coupled to a string via a massless 
spring. The long-dashed lines refer to the equilibrium positions of the string 
and the oscillator. The inset shows the three forces that cancel at the massless 
attachment point A {x = 0): the left and right tensions of the string and the 
elastic force of the coupling spring. 

respect to its equilibrium position will be denoted by X{t). 



2.2 Equations of motion 

The equation of motion of the oscillator is 



M- 



di2 



+ MnlX^~K{X~io) 



(la) 



where £,Q{t)^£,{0,t) is the displacement of A and fio = ^JK/M is the harmonic 
frequency of the free oscillator. For a; 7^ 0, ^ fulfills the one-dimensional d'Alem- 
bert equation 

5^e la^e . (^,) 



dx^ c^ dt^ 







where the wave velocity on the string is given by c=^/Tq/]jlq. The derivation 
of (fTb|) is a major step in every introductory course on waves (Crawford, 1968, 
§ 2.2, for instance). Less common is perhaps the refinement of sticking the string 
to a "mattress" (see figure [U made of n massless springs per unit length along 
the 2:-axis, each of them having a stiffness k. The restoring force per unit length 
due to the mattress, — n>ir^, turns the d'Alembert equation into a Klein-Gordon 
equation 



yg 1 d'S, 
dx"^ 






(x^O) 



(lb') 



where uJo^C\/i^^/To. This equation governs also the propagation of the elec- 
tromagnetic field in a rectangular waveguide (Feynman et al., 1970, chap. 24) 
and is also the relativistic equation of a quantum particle of mass Hujq/c^, c 
being the velocity of light in the vacuum. 




Figure 2: A mechanical model that leads to a Klein-Gordon equation: the 
string is attached to an elastic support that can be seen as a tight collection 
of massless springs uniformly distributed per unit length along the string. For 
another mechanical example, see (Crawford, 1968, § 3.5). 



Since the attachment point A is massless, the sum of the three forces ap- 
plied to it vanishes as depicted in the inset of figure [T] To first order in the 
slope d{x,t) — d^/dx{x,t), the transverse component of the right tension ap- 
plied to A is given by To^(0"'",i) the limit when x -^ keeping x > 0. On 
the other side, the transverse component of the left tension is — To'i?(0^,<). The 
restoring force of the coupling string corresponds to the opposite of the left hand 
side of (|lap . namely k{X—^q). Therefore the coupling introduces a discontinuity 
of the slope of the string at a; = 0: 



ro|i(0-,,) 



T.|,o-,<) 



-K(X(t) - &(i)) 



(Ic) 



2.3 Dimensionless quantities 



The fundamental units will be chosen to be T=M/{fj,oc) for the times, L = M/ jjlq 
for the lengths and M for the masses. The model is therefore uniquely deter- 
mined in terms of dimensionless quantities defined by the appropriate rescaling: 

Xeff=x/L, t,s"^t/T, Ccff^^/i, Xeff=X/i,^r!o^eff = f^o/(r-l), A^eff = ^/(Af T-^), 

etc. These conventions correspond to CcS^c/{LT^^) = 1. There is no need for 
any quantization as long as the effective Planck constant hcs = h/{ML'^T^^) — 
hljLo/{M'^c) ^ 1. For simplifying the notations, in the following we will drop 
the "effective" subscript and work directly with Af = 1, c= 1, /io = 1 and Tq = 1. 
Introducing the shifteco frequency 



n, ""^^ ^ni + K , 

the equations governing the dynamics of the system become 

{0-,t)^K{^o{t)-X{t)). 
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dx'^ 



and 
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(3a) 
(3b) 

(3c) 



^With more pedantry, one could speak of the ronormalized frequency. Here the shift can 
be understood by the effective restoring force —{K + k)X that M actually feels. 



with the help of the Dirac distribution S, equations (|3b[) and ([5c)) can be com- 
bined in one equation, valid for aU x, 

0-0--oe = '^'5(x)(^-X). (3bc) 

Indeed, ([5c)) is recovered after integrating ([3bcP between x = —e and x~+e 
when e ^ 0^, since ^ and its time derivatives are continuous (and finite) every- 
where. 

2.4 The Hamiltonian 

If one wants to prepare the ground for some perturbative treatment of some 
non-linear corrections, if one wants to quantize the model and/or to couple 
it to a thermal bath, one possible starting point is the Hamiltonian of the 
system expressed in term of some canonical variables. The continuous part of 
the system (the string) corresponds to a Hamiltonian density (Goldstein, 1980, 
§ 12-4) involving a pair of canonically conjugate fields (^n{x,t),^{x,t)) whereas 
the oscillator is described in terms of {Px,X). The Poisson bracket has also a 
mixed structure of continuous and discrete variables: For any two Oi, O2 that 
are functions of {Px,X) and functionals of (tt,^), 

.„ ^1 dcf dOid02 d02dOi , ffSOid02 S02S0i\, .^. 

The Hamiltonian of the whole system generates the evolution of any func- 
tion(al) O via dO/dt — dtO + \H, OJ. For our system we have 



H "^^ ip| + if7^x2 + i.(X-eo)^ + i| 






da; ; (5) 



the corresponding Hamilton equations yield directly to ([3]) . The interaction term 
is completely different from the one used in a recent one-dimensional pedagogical 
model (Boozer, 2007). The latter emphasizes the recoil effect of the field on the 
mass Af (unlike in the present work, the external (translational) degrees of 
freedom for M are considered) whereas we are more interested in the resonance 
effects. 

3 Scattering 

3.1 Free modes 

As for any vibrating system, a normal mode is defined to be a particular collec- 
tive motion where all the degrees of freedom oscillate with the same frequency. 
In the absence of coupling (k = 0), for a given frequency w ^ 0, one can choose 
two independent normal modes (the free modes) on the string given by 

C|,(x,i) = -^e'(±^-"-*). (6) 



The wave number is obtained from the dispersion relation of the Klein-Gordon 
equation: 

k{uj) =^ y^w2 _ cj2 ^=^ uj{k) = ^Cj2 + fc2 (7) 



When Lu < ujQ, k is purely imaginary with a positive imaginary part and does 
not correspond to any traveling wave. When lo > loq, k is real positive and the 
two modes are two monochromatic counter-propagating waves. 

A real field ^ obeying the Klein-Gordon equation (|3b|) carries an linear 
energy density pe — ^[{d^/dt)^ + {d£,/dx)^ + Wq^^] ^^'^ ^ linear density of 
energy current je = —dS,/dtd^/dx. The conservation of energy takes the 
form dpe/dt + dje/dx = for x 7^ 0. For a monochromatic traveling wave 
of complex amplitude a, i{x^t) = ae'*-^'^^""*-', an elementary calculation shows 
that the average current over one period is given by 

{je) = ±fcw|a| 2/2 (fc real), (8) 

whereas for an evanescent wave 

(je) =0 (fc imaginary). (9) 

3.2 Definitions of the in and out asymptotic modes 

When a coupling is present (k > 0), the two free waves ^ are not solutions 
of ^ anymore. A monochromatic traveling wave will be partially reflected 
(resp. transmitted) by the oscillator with a reflection (resp. transmission) 
coefficient p (resp. r) that is a complex function of k or oj. The linearity of 
equations ^ guarantees that the frequency will be unchanged by scattering. 
Indeed, to describe such a scattering process, a relevant choice for the two modes 
of frequency w is to look for in-states, defined (for real positive k) to be of the 
form (see figure [3] a)): 

1 ( i(kx-t^t) I i(-kx-ujt) for T < n ■ 

^^ ' v/2^ Ire'e^^-"*) forx^O; ^ ' 



and, since the scattering is symmetric with respect to a; i— > — a;, 

A f 1 fT-pi(-fea;-cjt) for T < n ■ 

with, for both modes, the same amplitude xi^) for the oscillator that can be 
interpreted, using the language of linear response theory, as a susceptibility: 

X-(i)=x(t^)e-''^*. (11) 

These modes form a set of waves that is suited for constructing localised wave- 
packets that look like free wave-packets when i — > —00 (i.e. far away from the 
oscillator). For instance, when considering a localised wave-packet traveling to 
the right, we shall havqj: 

cpik) e|(x, t) dk ^-— . J ifiik) a^ix, t)dk^^J ^(k) e'C^^-'^dA: . 

(12) 



■^This can be understood with the stationary phase approximation. If •f> is concentrated 
around fcp > 0, the wave-packets travel with the group velocity ±da;/dfe(feo) ^ 0. When t — ► 
— oo, the dominant contributions to the integrals in H12|l form one wave-packet located in x < 
and traveling to the right. 
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Figure 3: Several choices for the normal modes with frequency uj. The wavy 
arrows symbolise monochromatic traveling waves whose complex amplitude is 
specified below, a) the in-state given by (jlOap : b) the in-state given by (jlObp : 
c) the out-state given by (|15ap and d) shows the general coefficients that will 
be linked by the S ot T matrix defined respectively by (PTj) and 



The continuity of ^ at a; = implies that 

1 + p = T 



(13) 



Conservation of energy implies that the average energetic current is conserved 
in the stationary regime. From ([8]), we must have 



1 



\P\' 



(14) 



Moreover, the equations are real and invariant under the time reversal, therefore 
ii ^{x,t) is a solution, so is its complex conjugate (^(x,t))* and £,{x,—t). If we 
consider the solution (C™(2^j~^)) ; then we obtain the mode depicted in [3] c) 
obtained from[3]b) by reversing the orientation of the arrows and by conjugating 
the amplitudes. This procedure defines the out-modes that behave like free 
modes for remote future times when packed in localised superpositions. We will 
have (see figure [3] c)) 



a"*(^,o'= (en^,-<; 



2tt 



^* ^i(kx—Ljt) 
^i(kx—ujt) 



* J\{ — kx—ujt) 
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dot 



e^\x,t)- unx,-t) 



27r 



P e 



Ji{ — kx—ujt) I „* A{kx—Ljt) 
^* ^i{-kx-iAjt) 



for a; < ; 
for X ^ , 

(15a) 



for a; < ; 
for a; ^ ; 



(15b) 



with 

-ILOt 



X™*(t)=x*('^)e""^* . (16) 

If we interpret ^™' as a superposition of in-modes coming from both sides that 
conspire to product no wave traveling to the right for a; < 0, we get 

T*= ^; p* = ^. (17) 



More mathematically, £^°^*^{x,t) can be seen as the continuation of £,™{x,t) 
to the domain of negative fc's. Indeed, we have ^'^"*(2:,t) = ^""^,(2;,^) provided 
that we define 

ri-k)''^'{T{k)y; pi-k)"^' (pik))* . (18) 

3.3 Definitions of the scattering and transfer matrices 

The in-modes and the out-modes are two possible bases for describing a scat- 
tered wave-packet. These bases can be obtained one from each other by linear 
transformations; the linearity of the equations of our model implies that they 
connect waves with the same frequency only, which is a major simplification. 
A typical scattering experiment consists in preparing one wave-packet traveling 
towards the scatterer (the oscillator). Long before the diffusion, this ingoing 
wave-packet is a simple superposition of in-modes. Long after the diffusion, 
we get two outgoing wave-packets that are naturally described in term of out- 
modes. The passage from the in-basis to the out-basis is described in term of 
the scattering matrix 5* that encapsulates all the information about the possible 
scattering processes It is made of 2 x 2 blocks S'(w) defined by 

(19) 

The decomposition of each out-mode in term of the two in-modes for, say, x ^ 0, 
leads to 

P+T p + T / 

In other words, for the general monochromatic wave 

^(^'^) = |c,eH— )+i.,eK— ) forx^o! ^''^ 

the S matrix connects linearly the coefficients: 




^In the literature, specially within the context of scattering of quantum waves (Taylor, 1972, 
§ 2c, for instance), the matrices that connect the free waves to the in- waves on the one hand 
and the free waves to the out-waves on the other hand are often introduced under the name of 
M0ller operators with the caveat that unlike the free states, the set of scattering states may 
be incomplete, that is insufficient to construct all the states. As we will see in § 14.21 to get a 
complete basis one may add to the in-states lllOI I the bounded modes when existing. 



In the absence of scattering (r = 1, p = 0), 5 simply reduces to the identity. 
The unitarity of S, which can be checked on (PU]) . can be seen as a direct 
consequence of the conservation of energy since, from ([SJ, the norm of the two 
vectors involved in ^ is preserved: |C+p + jD-p = |C_p + [ZJ+p. 

If one wants to calculate the diffusion by several scatterers, it is more con- 
venient to introduce the transfer matrix T whose 2x2 blocks are defined to 
connect the left coefficients to the right coefficients, 

Then we have 

^H=r; ;j (24) 

whose determinant is one. The addition of one scatterer on the string corre- 
sponds to a multiplication by a T matrix. 

3.4 Physical interpretation of the solutions — Resonant 
scattering 



The definitions and the general properties presented in §§ 13.21 and 13.31 are valid 
for any non-dissipative punctual scatterer. As far as our model is concerned, 
inserting the expression ()10a|) with the oscillation (fTT|) in equations ^ yields 
to a linear system that can be solved straightforwardly: 



2.Ju^^^^ ^2 _ ^2 



p(a;)=-i + ie-2i')(-) = ^ , (26) 

2 2 ^_.V^f3^_^2_^2 



"o 



with 



and 



def 






rjiu;) - arctan | „ , , ., -^^ ) (27) 



^(^) = -2 — 2 -^ — n — iT • (2^) 

Even though the coefficients p and r were first defined for traveling waves, 
i.e. for Lu > loq, the above expressions can be continued for uj < loq. We will 
understand the physical interpretation of this procedure when we will study 
radiation in §|31 As long a.s ui > luq, ([H]) . pT)) and p^ hold. We can check that 
the ultra-violet limit w ^ cxd is equivalent to the limit of weak coupling where 
the oscillator becomes transparent: r — > 1 and p — *■ 0. Another case where the 
coupling is inefhcient is when A and X both oscillate in phase with uj — ^Iq 
since the coupling spring remains unstretched. More interesting is the resonant 




fin n^ 




Figure 4: Graphs of |t(u;)| (dashed Ime) and \p{uj)\ (sohd line) given by 
and (|26p when wq < ^o- One can observe a resonance scattering spike for uj — Qk 
and an antiresonance for a; = fi^ where no scattering occurs. The quahty factor 
is a) Q ~ 15 and b) Q ~ 200. Let us mention that when flo < ojq < i^K, the 
antiresonance has vanished and the local maximum |p(r2fe)| = 1 is too soft to 
be called a resonance. When fi^ < wq, no scattering resonance occurs and \p\ 
decreases monotonically from 1 at uj ~ ujq to when uj — s- +cx). 



scattering that occurs, provided that luq < ^k, when the ingoing wave that 
forces the oscillator has precisely the same frequency as the shifted frequency 
of the latter, that is cj = flk- Then, the scattering is the most efficient since no 
transmission occur (p = — 1, t = 0). The resonance spike can be seen in figure |4| 
and its quality factor can be evaluated from its width Atu when |/3{r2K±Aa;/2)| = 
l/\/2: for a small coupling. 



dcf 



n. 



2ni 



7e:-Ti^V"o--§(i + oW) 



Aoj 
and therefore, the smaller k, the better the quality of the resonance. 



(29) 



4 Radiation, damping and bounded mode 

The general idea that damping and therefore irreversibility emerge because of 
the interaction with a large number of degrees of freedom can be illustrated 
explicitly on our model. If we choose initial conditions such that the spring is 
at rest at i = 0, the entire energy being contained in the oscillator, for instance 



e(x,o) = o 



dt 



(x,0)=0; X{0)^Xa; X(0) = 



(30) 



the energy transfer to the string will damp the oscillations of M and the latter 
may completely lose its energy far before the energy can get back from the string 
if its boundary is far away from the oscillator (for a string of length i, Poincare 
recurrence time is of order i/c if there were no dispersion). The equation of 
motion of the oscillator is particularly simple when the string is non-dispersive 
(wq = 0) and therefore we will start by studying this case. However, we will also 
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consider the case of the Klcm-Gordon string because when ujq > fioi we will see 
that there exists a stable mode of the oscillator at a frequency w^ > whose 
dissipation is blocked because u>b < ujq. Its vibration does not decay because 
at this frequency, only evanescent waves can exist on the string, which do not 
carry away energy current on average (see ^)- 

4.1 The d'Alembert string 

With the initial conditions ([50]) , the general form of the radiated waves on the 
string will be ^{x,t) — (^o{t — \x\): each of the two wave-packets travels away 
from X ^ without distortion when luq = 0. Equation pep becomes 



2^0 + <o = ^X (31) 

and the elimination of ^o from (pij) and (Pa|) yields to 

on'^ . 2 
X + nlX = ^X--X. (32) 

K K 

The two terms in the right hand side are dissipative forces. The first one has the 
familiar taste of the viscous resistive force whereas the second has the flavour of 
the Schott term 2e^ 'i'/S (Rohrlich, 2000, eq. (2.7b)) in the Abraham-Lorentz- 
Dirac equation which governs the dynamics of an electric charge e that takes 
into account the electromagnetic self-force of the charge. The major difference 
is the sign of the coefficient in front of the third-derivative. Unlike the Schott 
term, the negative sign in ((32|) prevents the spurious exponentially accelerating 
solutions. It can be clearly seen that the irreversibility due to dissipation comes 
straightforwardly from the choice of initial conditions ((5(11) that break the time- 
reversal symmetry under which the original equations are invariant. 

Many models of an oscillator coupled to one-dimensional waves are recov- 
ered in the limit of strong coupling k —^ -|-oo (see the references given in the 
introduction, for instance when the mass is directly attached on the string). In 
that case, only the viscous force remains in (|32p and we immediately get the 
well-known damped oscillator X + 2X + HqX = 0. 

Looking for exponential solutions X{t) = X(0)e^* yields to the characteristic 
equation of ([32]) : 

z^ + -Kz^ + nlz + -Kfll = . (33) 

There are three solutions, one real zq and two complex 2;+,z_ all having a 
strictly negative real part. Perturbatively in k, we have 



2 2^2 



9 



For generic initial conditions, including (|30p . where the string is at rest the 
energy of the oscillator will exponentially decay like e~'"* at the rate 
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In the language of particle physics, the stable non-interacting particle (the mode 
of the free oscillator) has been destabilised into a metastable particle of life- 
time r~^ because of its interactions. 



4.2 The Klein-Gordon string 

When loq > 0, one cannot get a differential equation for X{t) but must keep 
working with its temporal Fourier transform 



X(cj) =^^ / X{t)e"^'dt, 



(36) 



together with a superposition of purely radiated waves of the form ^(x, t) 
(\/2^)"Vl('^)e'^''''''"'^*^dw. Inserting them in ^, X must satisfjH 



k{uj' - n'o) - 2iJuj^ - Luliu^ - ni) 



X{lu) = 



(37) 



Therefore X vanishes everywhere but at the frequencies that cancel the brackets. 
These are precisely the poles of r and therefore of p = t — 1 given by ([^5]) 
and (|26p . Indeed, for pure radiative modes, the ingoing waves vanish (C_ — 
_D_i_ = 0) and therefore the matrix element T22 must go to infinity in order to 
keep D_ finite (see figure [Sid) and equations ([25]) and ([M]) ). Letting Z = —uP'., 
we look for the solutions of the cubic equation 





(^ + 


ujI){Z^VlI 


Perturbatively 


in K, those are 






^0 = 


-^0 + 


4 




^+ = 


z*_ = - 


-^l 



+ nf-2^{Z + nl) 



nlf = . 



mi - ^i) 



2.M 



0(/.4) ; 



(38) 

(39a) 
(39b) 



When wq — ^ 0, we recover Zq — > Zg and Z± —> z^. The physical frequencies 
will be the three square roots i^/Zo,iy/Z± whose imaginary part is not positive: 
The typical decay rate of energy will be given by the nearest root a;,nin to the 
real axis: F = — 2Im(a;inin)- As long as 2ujo < k <^ fl^, all the three frequencies 
have strictly negative real part. The decay rate is given by 



F = 



rio 



Q 2noy/ni-, 



0(^3) 



(40) 



It is a very general feature that the poles of the S matrix are associated with re- 
sonances and, more precisely, that their imaginary part provide the decay rates, 
which are proportional to the inverse of the quality factor of the resonances. 

When K < 2ujo, Zq is negative, one residual oscillation persists at fre- 
quency ajf, = \/— Zq. For uji, < ujq, no transfer of energy is allowed; only evanes- 
cent waves are created and those do not carry any average energy current. Un- 
like the scattering states, this non-decaying mode is spatially localised. More 



*The presence of the square root in l|37p is tlie reason that prevents us from obtaining a 
local differential operator for X{t). 
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generally, any bounded mode has a purely real frequency ut that must be less 
than luq since Zf, + uJq = uiq — uj^ > in order to fulfil ((55)) . From ([7]), fc(wfc) 
is therefore purely imaginary. Moreover, in order to cancel the bracket in p7p . 
ujf must lie in between Qq and il^. For simplicity, let us introduce the auxiliary 
parameter T = (riQ — (jJq)/k and the real positive variable u^\k\/^/K; a stable 
mode will exist if the cubic equation 

2uiu^ + T + 1) + ^.{u^ + T) = . (41) 

has a positive real solution. This can be achieved for T < only i.e. in a 
regime where ilo < wq. For T < 0, the product of the roots of the left hand side 
of (liTj) . u+U-Ub is —\/kT > 0. If two roots are complex conjugated, therefore 
the third one is necessarily positive. If the three roots are real, either only one 
is positive or all three of them are. The latter case must be ruled out since the 
sum u+ + w_ + Uh = —^/k is strictly negative. We have therefore proved that 
a sufficient and necessary condition for a stable mode to exist is that wq > f^o. 
Its frequency is given by 



ujb = \/^o ~ ^""b (42) 



where Ub is the unique positive real solution of ()41|) . Perturbatively in k, we 
have 



^^ = »o + TT^ - 1° .^., ° J, ^ + Q(^^) (43) 



and the corresponding bounded mode is given by 



^\x, t) = Cb e-V'^o-'^'l^le-''^"* ; (44a) 

Xb{t)^j^^e-'-^r (44b) 



ni - ul 



The choice of the normalization, 



-1/2 



g^=( r^^ + (r,2"\..2 ] (45) 

will be justified below (equation [54]) . 

5 Some like it diagonal 

5.1 Normal coordinates 

What makes the model completely tractable is of course that it remains linear. 
However, the direct diagonalization of the quadratic Hamiltonian ([5]) remains 
particularly difficult. In that case, the trick is to solve the equations of motion 
to determine the normal modes first — this is precisely what we have done in 
the previous paragraphs — and then write the Hamiltonian in its diagonal form. 



H 



9 51 (^" + '^a'?^) ^^^o^al a^ (46) 
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in term of some (real) canonical coordinate {pa, (la}a or (complex) normal co- 
ordinates {a„}a that are associated with modes labelled by the discrete and/or 
continuous index co. We have 

^ ' ; (47) 




Qa = ^=={a*a + aj ; Pa ^ U ^{a^ - aj ; (48) 



and, for each pair {ai, a2}, 

{aai,aa^} = 0; {0^,1,0*^} = i^ai.aa ; (49) 

\Pai,Pa2ff — 0; \Qai,Qa2l — 0; \Pai,Qa2l — "qi,Q2 i (^O) 

where (5 stands for the Kronecker symbol or the Dirac distribution. The second 
step consists in determining the canonical transformation that expresses Oq in 
terms of some a priori known normal coordinates, namely some free normal 
modes aft-^Q. In our case this transformation is linear and will be transposed 
directly into the quantum theory by replacing the complex number a^ (resp. 
a* ) by the creation (resp. annihilation) operator a^ (resp. its Hermitian con- 
jugate a* ) of the ath one-particle eigenstate whose energy is fkOa- As we have 
seen, all the scattering states are twice degenerate, in the sense that each nor- 
mal frequency uj is associated with two independent states labelled by k and —k. 
These modes both diagonalize the Hamiltonian ([5|) . An infinite number of pairs 
of eigenvectors can be chosen to constitute a basis, among them, the in and 
out-states, which are particularly relevant as soon as we get into a quantum 
field theorjO. But it order to get ([H]) properly one must check that the set of 
modes is actually complete — i. e. that any kind of motion of our system can be 
described as a linear superposition of modes — and orthonormalized correctly 
in order to deal with canonical complex coordinates. Fourier analysis assures 
that the free states ^ constitute a complete set for describing the waves on the 
string. When interacting with the oscillator, if ojq > ^o one bounded state ex- 
ists that must be added to the in-modes (or to the out-modes) to get a genuine 
basis. Then, including the normal coordinates At of the bounded mode if there 
is any, (|46|l reads 



H ^ ujbAlA^ +J Lj{k) al,{k) a;^{k) dfc = ^6 A^A + / Lo{k) al^,{k) a^^,{k) dfc . 

(51) 
We chose the convention that, when fc > 0, ain(fc) (resp. aout(fc)) is constructed 
from ^'" (resp. ^™*) while ain(— fc) (resp. aout(— fc)) is constructed from ^™ 
(resp. er). 



^To avoid ambiguities we will often subscript the brace describing a set like {. . . }a^A to 
recall which indices are running and what is their range A if the latter does matter. 

® The normal coordinates ain(fc) and aout(fc) constructed from the scattering modes, once 
quantized into a\^{k) and aout{fc), allow the interpretation of the quantum states in term of 
asymptotic (quasi-)particles; more precisely the linear transformations from the free afr(fc) 
provide the explicit connection between the non-interacting states (the Fock space for bare 
particles including the free vacuum) and the interacting states (the Fock space for dressed 
particles including the interacting vacuum). 
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5.2 Orthonormalization 

Having a complete set of modes does not guarantee that they are orthogonal. 
Indeed, it may happen that two eigenvectors having a common eigenfrequency 
are not. For instance, one must check in one way or another that the modes (fTO|) 
are orthogonal and properly normalized. If we denote by S(i) a classical state 
represented by the displacement X{t) of the oscillator and the wave ^(x,i) on 
the string, the scalar product between two states Si(i) and S2(i) is 

Si(t).S2W=Xi*(t)X2(i)+ [ Ci{^,t)i^(x,t)dx. (52) 



It is shown in the appendix that if ^jj" (resp. 'B}\ ) stands for the mode S,^)^{x, t) 
(resp. £,'l^{x,t)) both with X™{t) ~ x(a;)e~"^*, then we have, for any (positive 
and/or negative) real pair (fci, fc2), 

Efjt)-E^{t) = d{k,-k2). (53) 

It is easy to see that we chose the normalization (|^5|) in order to get 

E\t) ■ E'it) = 1 . (54) 

The free states s![^. represented by X = and ([6]) are clearly orthonormal- 
ized, S|' • S|' = S{ki — k2), and form a complete basis if we add the state 
that allow to describe the motion of the oscillator, namely S^^^, represented 
by X = 1 and ^ = 0. 

5.3 The real symmetric modes 

The potential in ([5|) is a real definite positive symmetric quadratic form and 
therefore can be diagonalized in an orthogonal basis of real vectors. The natu- 
ral choice of retaining the real or the imaginary part of SlJ}j, actually provides two 
real modes but that are not orthogonal. A way to assure that we deal with an 
orthogonal basis, is to pick up a symmetry, say the parity x i— > —x, of the Hamil- 
tonian and classify the eigenmodes accordingly. The bounded state, if there is 
any, remains even. The in and out modes are not symmetric under space inver- 
sion but it is straightforward to obtain eigenmodes that are also eigenvectors of 
parity. For any complex factors c±, the combinations c±(S)5' ± S'"^,) are sym- 
metric/antisymmetric eigenvectors at any time with the eigenvalues ±1. After 
some algebraic manipulations using the expressions (|26p of p in term of rj given 
by (|27p . the symmetric and antisymmetric modes are represented, for /c > 0, by 

c+{a\x,t)+ei:(x,t)) = c+c-'^\^cos{k\x\-v)e-'^'; (55a) 



-{m^,t)-^fix,t)) = ic_y|sin(fcx)e^ 



(55b) 



The choice c+ — e"'/\/2 and c_ = — i/-\/2 leads to the real normalized symmetric 
modes, defined for fc > by 

'Z^^ = R{cj) (^rlj . (56) 
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with the unitary matrix 



/ e''''"-' e"''-"-' 



Constructing the real (anti)synimetric states from the out-modes at any time 
leads to the same Sj, since the common eigenspace to H and to the parity is of 
dimension one. To sum up, for fc > 0, S^^ is represented by 



1 r.,| „/, ,M„-iu.t. v+U\^ « C0S[77(W)] i^t 



it{x, t) = ^ cos[k\x\ - ri{Lo)] e-'-* ; X+{t) = __i^ e""^* ; (58a) 

^^ (a;, t) = ^ sin(fca:) e^''^* ; X" (t) = , (58b) 



and {Sb(i)} U {S^. (i)}fc>o is, at any time, an orthonormalized eigenbasis of 
symmetric or antisymmetric eigenvectors of the Hamiltonian with eigenvalues 
given respectively by dH]) and ([7]). 

5.4 An atom in a closed cavity 

The explicit canonical linear transformation that connects the free canonical 
variables to the interacting ones (the in or out modes via the real symmetric 
ones) is beyond the scope of this article and will be given and extensively inter- 
preted in a future paper where we will quantize our model. As explained above 
(see § 12.41 and also the note[H]), this is really interesting and beyond a purely 
academic exercise only if one wants to switch to quantum theory and/or sta- 
tistical physics. The quantum linear transformation between ain(fc) and afr(fc) 
appears to be a generalised Bogoliubov transformation and our model provides 
an explicit construction of quasi-particles in terms of free particles. 

However, the real symmetric modes that have been founded in the previous 
section remain interesting at the less advanced level of the present article because 
they are the natural modes to work with when the finite size (. of the string 
becomes relevant. Indeed, when £/c is not too large compared to the typical 
time r~^ characterising the radiations of the oscillator, the discrete character 
of the spectrum of the non-interacting string can be "feeled" by the oscillator. 
When boundary conditions are imposed, say £,{£/2,t) = £,{—£/2,t) = 0, the 
discrete (even) spectrum of the whole system is modified by the presence of the 
oscillator and, from (|58a|) given by the {fc„}„gz that fulfill the equations 

1 TT 

- k^e - ■q{u}{kn)) = - +mT -^=^ tan {T]{k)) = tan(M/2 - 7r/2) (59) 

that can be solved graphically (figure [5]). The frequency Hq of the free oscilla- 
tions of the mass inserts in the spectrum of the string. The even spectrum will 
differ from the non-interacting case when c"* is significantly different from one. 
For resonances with high quality, it will not affect the frequencies that are away 
from the resonant frequency. 

What one gets here, for wq = is an elementary model of an atom in a 
(perfect) electrodynamics cavity of size £ (some imperfections can be taken into 
account if we relax the Dirichlet boundary conditions and put partially reflec- 
tives "mirrors" on the string). The field may or may not be quantized and. 
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^0 



O. 



Figure 5: Graphical resolution of equations ([591) that provide the even frequency 
spectrum when the oscillator is attached to the middle of string of finite size £. 
The frequencies, represented by small up triangles on an horizontal line at the 
bottom of the figure, are centered on the abscissae of the intersections of the 
graph lli i—s- tan [ri{u!)) (thick solid line) with the graphs u i— > tan {k{uj)i — 7r/2) 
(thin solid lines) for lo > loi^. As in figure[4]a), the quality factor of the resonance 
is Q ~ 15. The down triangles indicate the spectrum of vibration of the Klein- 
Gordon string alone of finite size I obtained for tan {k{uj)i/2 — 7r/2) — 0, that 

lY /P with n a positive integer. 



is for ujn — \A 



H2n- 



not to speak of lasers, we obtain a sort of primer for the widespread physics of 
quantum electrodynamics cavities that have been realised in laboratory to test 
successfully some fundamental concepts in quantum physics (Haroche & Ray- 
mond, 2006). The purely mechanical model for infinite k (the mass is directly 
attached on the string) has been carefully studied with experiments in (Gomez 
et al., 2007). 

6 Conclusion 

In addition to a more detailed study of the finite size effects, another natu- 
ral development of the present work would be to deal with multiple scatterers. 
For instance, when there are two identical scatterers with wq > f^o, we expect 
that the degeneracies of the two bounded modes is broken and that the split- 
ting between the symmetric and the antisymmetric bounded modes decreases 
exponentially with the separation of the oscillators. Starting with initial condi- 
tions where only one oscillator has some energy, the beating between the two 
oscillators is an example of tunnelling due to the presence of evanescent waves 
connecting the two oscillators. 

Even before we quantize the whole system, our model may be interesting to 
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keep the field classical whereas only the oscillator is quantized. It would provide 
an illustration of say, the Fermi golden rule within the context of time-dependent 
perturbation theory (Cohen- Tannoudji et al., 1980, chap. XIII). However, we 
have proven that this golden rule transpires in our classical model since the 
transition rate (j35|l to the continuum of the modes is proportionnal to the 
square of the coupling strength k to first order in the perturbation. 

As it has been demonstrated, this model captures many fondamental phe- 
nomena that are important in many areas of physics and offers wide possibilities 
for pedagogical use. Above all, I hope it will help the reader to discover and/or 
to understand them more deeply. 
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7 Appendix: Normalization of the modes 



The construction of the real symmetric modes presented in § 15.31 leads to an 
orthogonal basis {Eb{t)} U {S^ {t)}k>o- When it exists (wq > ^o), the bounded 
state Sb has a norm unity. That they are eigenvectors for different eigenvalues 
of H or of the parity guarantees that S^ • S^.^ oc 6{ki — ^2) and this appendix 
proves that the proportionality factor is indeed unity. 

Rewriting (jlOap with the help of p3p and with the Heaviside step function 8, 

C^{x, t) - -^ [e'^-^ + p{k) e-'^-Q{-x) + p{k) c^^-Q{x)\ e""^* , (60) 

we have, with pi = p{ki) and p2 — ^(^2), 



/ 



CkA^,t) ^kA^,t) dx ^ 6{ki-k2)+ ' 



2tt 



P2 



1 1 



1 1 \ , MP2 

Pi 



k2 + ki+ iO+ k2 - ki+iO 

(61) 



^ fc2 - fci + iO+ k2 + ki - iO+ / fc2 - fci + iO 
We have used the identity (0+ stands for the limit e ^ keeping e > 0) 



+00 _ jgifexo 



^'"^^=M:ioT (62) 



iTrS{k) (63) 



ixo 

valid for any real k. The other identity 

k + iQ+ k 
allows to convert (|6ip in terms of the 6 distribution and of the Cauchy principal 
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value p: 

Cfci (2;, t) 4^ (x, t) dx = S{ki ~ /C2) 
(P2-PD7- 



+ pt'^(fci+fc2) + — 

ZTT 



P 



fcl 



+ (pI +P2+2ptP2)- 



P 



fcl 



(64) 



In fact, the coefficients of tfie principal values both vanish when the respective 
denominators cancel (we use psp . ([T^ and (fH|) . then p + p* + 2\p\'^ = follows) 
and we can drop the symbol p. The S{ki + fc2) can also be forgotten for, to 
constitute the basis, we retain only strictly positive values of ki and k2. A little 
bit of algebraic jugglery with (PS)1 and ^E^ allows to check that 



1 

2^ 



P2- Pi , P*l+P2 + 2plp2 



h 



ki 



2^{ni-ui){ni 



77 = -X*{^l)x{^2) 

(65) 

with T„ = T{kn), oJn — ^(fcn) {fi — 1, 2). Then we have proved that, for fci > 
and ko > 0, 



4" it) ■ 4" it) = (^i" (0) * ^s (t) + / a, i^, t) 4. {x. t) dx = j(fc 



(66) 



The space inversion of this identity leads immediately to S"^. it) ■ :i™^ [t) = 
S{ki — k2). At last, with the same techniques we can obtain 



Ck, i^, t) 4^ (x, t) dx = (1 + pi)S{ki + fcz) + 



1 
2^ 



2^2 Ti* -T2- Ti* 
/C2 — ki 



fc2 + fci 



(67) 



As above, for ki and ^2 both strictly positive, 5{ki + ^2) vanishes whereas the 
second term on the right hand side is precisely — (X"(i)) X]^f. (t). Therefore 

E^ (t) ■ SL\ it) = {Xj^ it)) * XX it) + I Ck, (^, t) 4. ix,t)dx^O. (68) 



The complex conjugation and the time reversal t 1-^ —t oi the above relations 
allows to show that the out-modes are also orthonormalised. The orthonormal- 
ization of {S^, }fc>o follows from (fSS)) and ((57)) . To sum up, we have obtained 
the following scalar products, for any fci > and ^2 > 0: 



S±\W-S'i\,(i)=<5(fci-fc2) 
SS\(i).SSUi)-'5(fci-fc2) 



ss\(i)-s™4(i) = o 

Elit)-Elit) = 0, 



2^,W-S±(i) = <5(fc2-fc2) 
and, when loq > flo for a unique bounded state to exist, 

Ebit) ■ Etit) = 1 ; 

S!t\^ it) ■E,{t)=Q; ElX it) • S, (t) = ; S± (i) • 5, (t) 







(69) 
(70) 
(71) 

(72) 
(73) 



Three eigenbases for the Hamiltonian have been chosen: {Sfc(i)} U {^)f (t)}fegR, 
{Ebit)} U {S^^*(i)}feeR and {S,,(i)} U {S^(i)}fc>o. The passage from one to the 
other is done with unitary matrices made of independent 2x2 unitary blocks 
of S'(aj) or i?(aj) given by (|^ and ([57)1 respectively. 
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